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Abstract 

A central feature in the Hilbert space formulation of classical me- 
chanics is the quantisation of classical Liouville densities, leading to what 
may be termed term Groenewold operators. We investigate the spectra 
of the Groenewold operators that correspond to Gaussian and to certain 
uniform Liouville densities. We show that when the classical coordinate- 
momentum uncertainty product falls below Heisenberg's limit, the Groe- 
newold operators in the Gaussian case develop negative eigenvalues and 
eigenvalues larger than 1. However, in the uniform case, negative eigen- 
values are shown to persist for arbitrarily large values of the classical 
uncertainty product. 

1 Introduction 

The phase space formulation of quantum mechanics has its origins in the work 
of Weyl [P and Wigner [5j more than 70 years ago. In this formulation, quan- 
tum observables are represented by phase space functions and quantum states 
by quasi-density (Wigner) functions that evolve in time according to a deforma- 
tion of the classical Liouville equation. Although many calculations are more 
difficult to perform in the phase space context, this formulation has the distinct 
advantage of the common language of phase space for the purpose of comparing 
classical and quantum dynamics. Indeed, it was for such a purpose (the calcu- 
lation of quantum corrections to thermodynamic quantities) that Wigner first 
introduced his now famous function. 

This begs the following question: if quantum mechanics can be described 
in the classical setting of phase space, can classical mechanics be described in 
the quantum setting of Hilbert space? In fact, it was already recognised by 
Groenewold [3] nearly 60 years ago that this is possible and, after related work 
in the interim [SJ EJ |Sj , it has recently been made explicit . The key to 
understanding this is to recognise the unitarity, and hence the invertibility, of 
the Weyl- Wigner transform |10l lllj , which maps quantum operators into phase 
space functions. 
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The inverse transform is Weyl's quantisation map, which transforms classi- 
cal functions (observables) into linear operators on Hilbert space. Just as the 
quantum density operator is mapped into a quasi-probability function on phase 
space (a Wigncr function) by the Weyl-Wigncr transform, so a classical Liou- 
villc density function on phase space is mapped by the inverse transform into a 
quasi-density operator on Hilbert space (which we suggest to call a Groencwold 
operator). Just as a Wigncr function is normalised on all of phase space, but 
is not in general a non-negative definite function, so a Groenewold operator 
has unit trace, but is not in general a non-negative definite operator. Just as 
quantum averages can be calculated in the phase space formulation of quantum 
mechanics using the Wigner function in the usual way of a density function, so 
classical averages can be calculated in the Hilbert space formulation of classical 
mechanics using the Groenewold operator in the usual way of a density opera- 
tor. Finally, just as the Wigner function evolves in time in accordance with a 
deformation of classical mechanical evolution, governed by a deformed Poisson 
bracket (the star or Moyal bracket), so the Groenewold operator evolves accord- 
ing to a deformation of quantum mechanical evolution, governed by a deformed 
commutator bracket 0, which has been called elsewhere the odot bracket jUJ- 

The time evolution of Groenewold operators in Hilbert space can also be 
described using a superoperator formalism that is particularly convenient for 
numerical work, as shown by Muga et. al. 0IH]. More recently [3], an expres- 
sion has been derived for the classical evolution on Hilbert space as a series in 
increasing powers of H with the quantum evolution as leading term. This form 
of the evolution suggests the possibility of probing the classical/quantum inter- 
face from the quantum side, by calculating classical corrections to a quantum 
evolution, rather than quantum corrections to a classical evolution. 

It is desirable to understand the spectra of Groenewold quasi-density opera- 
tors because of their primary role as the Hilbert space representatives of classi- 
cal states. Some insight into how negative eigenvalues of Groencwold operators 
can arise is provided by the recent observation [15] that if a nonlinear classical 
evolution is applied to an initial quantum density operator corresponding to 
a coherent state (such an operator is also a non-negative definite Groenewold 
operator, corresponding to a Gaussian initial density on phase space), then the 
evolved operator develops negative eigenvalues. This is a special case of a result 
obtained earlier [7] for initial density operators corresponding to more general 
pure states. With regard to characterising those Groenewold operators that are 
non-negative definite, we note that every non-negative Wigner function can also 
be regarded as a Liouville density, and the corresponding true density operator 
is therefore also a Groenewold operator that is non-negative definite. However, 
the problem of identifying all non-negative Wigner functions has been solved 
only partially |12| . The only pure states having non- negative Wigner functions 
are coherent states (modulo a linear canonical transformation), and all such 
Wigner functions are Gaussians. Although progress has been made in extend- 
ing these results to mixed states ^J^], the extent of the set of all non-negative 
mixed-state Wigner functions is still unknown. 

In what follows, we are concerned with properties of Groenewold operators 
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at a fixed time, and not with their evolution. We consider a classical system 
with one linear degree of freedom, in the phase plane T and in Hilbert space Ti. 
We examine the eigenvalue spectra of Groenewold operators in Ti corresponding 
to (A) classical densities that are Gaussian on T and (B) classical densities that 
are uniform on a circular or elliptical disk in T, and zero elsewhere. In both 
cases the eigenvalues can be calculated exactly. In Case (A), we show that 
negative eigenvalues, and eigenvalues that exceed 1, appear when the classical 
uncertainty product AqAp falls below H/2, the minimum allowed by quantum 
mechanics, and we discuss the limiting form of the spectrum as the Gaussian 
approaches a delta-function. 

The reader might be misled into thinking that a general classical density will 
give rise to a positive Groenewold operator whenever the uncertainty product is 
greater than H/2. This is not so, and we show that in Case (B) the corresponding 
Groenewold operators have negative eigenvalues for all values of AqAp. 



2 Groenewold quasi-density operators 

We denote a general classical (Liouville) density on T by p(q,p). It satisfies 

p(q,p)>0, ^p(q,p) = l. (1) 

The corresponding Groenewold operator may be defined through the inverse 
Weyl-Wigner transform: 

p = 27ThW- 1 (p) = J^p(q,p)A(q,p)dqdp, (2) 

where A(q,p) is the Weyl-Wigner kernel ^5], which may be defined as ^7] 

A(q,p) = 2D(q,p)PD(q, P y = 2A(2q,2p)P . (3) 

Here P denotes the parity operator and D(q,p) = exp(i(c[p — pq) / H) is a unitary 
displacement operator. Note that the eigenvalues of A(q,p) are ±2 an d 
since Groenewold operators can be expressed as a convex combination of these 
operators as in [21 their eigenvalues lie in the range [—2, 2], whereas those for a 
true density operator are restricted to the interval [0, 1]. 

If A(q,p), B(q,p) are classical functions on phase space and A = W~ 1 (A), 
B = W _1 (J3) are their operator images under the inverse Weyl-Wigner trans- 
form, then we have 

^ j^A{q,p)dqdp = Tr(i) , ^J^A(q,p)B(q,p)dqdp^Tr(AB). (4) 

Note also that if A(q,p) is real- valued, A is Hermitian. These properties ensure 
that a Groenewold operator p = W~ 1 (p) has a number of important features 
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in common with true density operators: it is an Hermitian operator with unit 
trace [5] and the trace of the product of two Groenewold operators satisfies 

Tr(pp') = 27rhJ^p{q, P )p'(q,p)dqdp > 0. (5) 

Furthermore, as mentioned above, the calculation of classical averages on Hilbert 
space takes the familiar form used for quantum averages, with the Groenewold 
operator acting in the role of a density operator: 

(A) = J p{q,p)A(q,p)dqdp = Tv(pA) . (6) 



3 Gaussian densities 



Despite the similarities between Groenewold operators and true density opera- 
tors, there are important differences, the most obvious of which is the failure of 
the Groenewold operator to be non-negative definite in general. 

The uncertainty product AqAp associated with a classical density is non- 
negative but can be arbitrarily small, and it is not surprising to find signatures 
of 'non-quantum' behaviour in the Hilbert space representation of classical den- 
sities appearing when this uncertainty product is less than ft/2, the minimum 
value allowed in quantum mechanics. 

In order to investigate this further, we introduce the class of Gaussian den- 
sities 

where /?, 7 are constants with units of Length and Momentum respectively. The 
corresponding classical uncertainty is given by AqAp = [3^/2. 

From (J2J we see that the corresponding Groenewold operators are given by 

Pfl = I P$(q, P )A(q,p)dr. (8) 
•/ r 

We now introduce the complex-conjugate variables a, a with a = (\J j/ j3 q + 



i\[Wpyp)l V2h and the corresponding annihilation-creation operator pair a,aJ 
with a = (1/7/ ' f3q + iyj l3/~ip)/\/2h, together with the "number operator" N = 
a^a. This allows us to define the Fock basis 

{\n),n = 0,1,2,...}, a|0) = 0, |n) = -^|0) , N\n) = n\n) . (9) 



(G) 

The reason for introducing this basis is that the Groenewold operators p^ ' 

commute with N and are diagonal on the Fock states. This can be seen by a 
direct calculation of their matrix elements: 



n\p^\m) = (m\pfa\ri) = / p^(q,p){n\A(q,p)\m)dqdp , (10) 



r 
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where, in terms of the variables a, a, (n|A(a, a)\m) is given by |18| 



(n|A(a,a)|m) = 2(-l) m W — (2a)"- m L"- m (4|a| 2 )e- 2|0!| , n>m. (11) 



If wc introduce the polar-type variables (t, </)) with a = \Jt/2cxp(i(j)) and 
a = s/t/2 exp(— i0), then we may express (111) II as 



— z(n— m) 



> e 



/rrjl f '■• ■ n v 

/— / dt(V2t)(«— )/ 2 L^- m (2t)e- (1+ 7/3 )t . 
• 7 2tt/37 ,/ V n! J 

(12) 

Except in the case that n = m, the integral over (j) vanishes, and hence, p^ ^ is 
diagonal on the Fock states, with eigenvalues 

/ i \n roc ft 

Ai G) (A7) = / dte-^W'L^t), n = 0, 1,2, ... . (13) 

Note that this result generalises in the sense that if / = W _1 (/), where f{q,p) = 
f{q 2 /P 2 +P 2 /j 2 )i then [f,N] = and the eigenvalues of / satisfy a formula 
that generalises l|T5)l. 

After evaluating [2"U| , we get the eigenvalues of in the simple form 

W(fin) = ^(^)\ n = 0,l,2,.... (14) 



Since AgAp = /?7/2 and no quantum state exists in which AqAp < K/2, one 
expects that non-quantum features of p y p ^ should be observed for p7 < h. 
The result (|14|1 shows that this is indeed the case, since for (fa/ < H and n 
odd, the eigenvalue A„ 7) < 0, while at least one of the eigenvalues lies 
above 1. Hence for fij < H, when the classical density Pp^(q,p) may be called 

"strongly classical", the Groenewold operator py^ displays the non-quantum 
feature of negative eigenvalues and eigenvalues that exceed 1 . This is analogous 
to the negative values exhibited by Wigner functions corresponding to "strongly 
quantum states" , such as excited Fock states [21] ■ 

(G) 

When the operator p^ ^ is considered on the coordinate representation of 
Hilbert space, it appears as an integral operator with kernel |11| 

/oo 
pjg([s + y]/2,p) j^-^dp 
-OO 

= _J_ e -(*+y) 2 /(4/3 2 ) e -7 2 (*-y) 2 /(« 2 ) _ ( 15 ) 

(G) (O 

Because the An (/?, 7) are the eigenvalues of ph * , it follows that 



P^(x,y)ip n (y)dy = X^((3,^ n (x) , (16) 
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where f n {x) is the coordinate representative of the state \n). This has the 
well-known form of an oscillator eigenfunction, 

<p n (x) = const.H n { Xs f^/m) e^/WV , (17) 
where H n is the Hermitc polynomial [35], and then l|15f> and (|16fl give 

/oo 
e ^ +h)m . h)x . {0 , +H)yfm ^ Hn {y^m) dy = PV^^HP, 7) H n ( X ^Jph) . 
-oo 

(18) 

This is a known identity [201 an( i provides a check on the validity of (1141) . 

The spectral bounds on , are graphed in fig. ^ against the classical 
uncertainty AqAp (in units of n). When AqAp < h/2, we obtain the strongly 
classical Grocnewold operators discussed above. At the point AqAp = h/2, the 
Groenewold operator has only one non-zero eigenvalue Ao = 1, and is precisely 
the pure state density |0)(0|. For AqAp above this critical value, the spectrum 
is uniformly positive and bounded by below and by A = 2h/ (2h + [3-f) < 1 

~( G) 

above. As such, for (3^ > h, the Groenewold operator ^ may be considered 
as the mixed state quantum density 

oo oo 
n=0 n=0 



4 Uniform Densities 

In the Gaussian Case (A) above, there seems to be an obvious connection be- 
tween the value of the uncertainty product associated with the Liouville density 
and the appearance of negative eigenvalues of the corresponding Groenewold 
operator. This is not necessarily the case for more general densities, however, 
and we illustrate this point by considering Case (B), the family of uniformly 
distributed Liouville densities, defined on the interior of an ellipse in T, centred 
at the origin: 

f 4-, o< £ + 4<i, 

f TT \ 7T/37 ' — p 2 ~t Z — ' 

P%]{1,V) = (20) 
[ otherwise . 

The uncertainty product associated with p^(q,p) is AqAp = (3j/4. 
The corresponding Grocnewold operators are given by 

P%] = J r P { ^(q,p)A(q,p)dqdp. (21) 

We recognise these as scaled versions of operators that have been considered 
previously in a different context [22] - They also commute with N as defined 
above and their eigenvalues may be expressed as 

X "^ = P^-J rt «(2i)*, n = 0,l,2,.... (22) 



G 




Figure 1: Spectral bounds on ^ (left) and p^ ^ (right). Note that when the 
classical uncertainty product exceeds h/2, the Groenewold operators for the 
Gaussians considered in Case (A) are non-negative definite. On the right we 
show that no such transition occurs for the Groenewold operators corresponding 
to the uniform distributions considered in Case (B). 



In fig. n the bounds on the spectrum of p^ ^ are graphed against f3j (in units 
of h), from which one sees that although the bounds decrease in magnitude as 
AqAp rises above h/2, the lower bound remains below zero. 

There is a simple physical explanation for this phenomenon: by constructing 
a uniform distribution over a finite subregion of V, one is restricting both the 
position and momentum of a particle to a finite interval. This is impossible for 
a quantum system because the configuration and momentum space states are 
related by the Fourier transform and it is known that the Fourier transform of 
a function with compact support is entire |24j . This argument applies not only 
to Case (B) but to any Liouville density with compact support in T: for any 
such density, the corresponding Groenewold operator cannot be a true quantum 
density and it follows that since it has a trace equal to 1, it must have at least 
one negative eigenvalue in its spectrum. However, we note that for the example 
considered here, the negative values decrease in amplitude as the area of the 
ellipse increases. 



5 Conclusion 

Just as Wigner functions play a central role in the phase space formulation of 
quantum mechanics, so Groenewold operators play a central role in the Hilbert 
space formulation of classical mechanics, which offers the possibility of new in- 
sights into the classical-quantum interface. It is therefore important to elucidate 
the properties of Groenewold operators. In the present work, we have consid- 
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ercd exactly solvable examples to illustrate some specific properties. We have 
seen that "violation of the uncertainty principle" and restriction to a compact 
subregion of T can lead to eigenvalues outside of the range [0, 1], but this can 
occur due to the influence of other factors. The emergence of negative eigen- 
values when Groencwold operators that are initially true density operators are 
allowed to evolve classically rather than quantum mechanically is a prime 

example. There is evidence QSl to suggest that when the effects of the classical 
and quantum evolution diverge rapidly (such as for a delta- kicked rotator), the 
magnitude of the negative eigenvalues is much greater than for slowly diverging 
evolutions (such as for the Duffing oscillator). In recent work, we have shown 
that large negative eigenvalues arise in the classical evolution of an initial co- 
herent state density [55] in an anharmonic potential [2H| that causes the phase 
space density to develop "whorls" . 
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